Nomenclature
1 1 1 , Inlet of the stator
2 2 2 , Exit of the stator
3 3 3 , Inlet of the rotor
4 4 4 , Exit of the rotor
Thermodynamic variables
Working fluid name
p 01 p_{01} p 01 , stagnation pressure at turbine inlet
h 01 h_{01} h 01 , stagnation enthalpy at turbine inlet
p 4 p_{4} p 4 , static pressure at turbine outlet
m ˙ \dot{m} m ˙ , mass flow rate
Turbine design parameters
ν = u 3 / v 0 \nu= u_3 / v_0 ν = u 3 / v 0 , blade velocity ratio
R = ( h 3 − h 4 ) / ( h 1 − h 4 ) R=(h_3-h_4)/(h_1 - h_4) R = ( h 3 − h 4 ) / ( h 1 − h 4 ) , degree of reaction
α 1 \alpha_1 α 1 , absolute flow angle at the inlet of the stator
α 2 \alpha_2 α 2 , absolute flow angle at the exit of the stator
r 1 / r 2 r_1/r_2 r 1 / r 2 , radius ratio across stator
r 2 / r 3 r_2/r_3 r 2 / r 3 , radius ratio across interspace
r 3 / r 4 r_3/r_4 r 3 / r 4 , radius ratio across rotor
A R S = 1 2 ( H 1 + H 2 ) / ( r 2 − r 1 ) AR_{\mathrm{S}}=\tfrac{1}{2}(H_1 + H_2) / (r_2- r_1) A R S = 2 1 ( H 1 + H 2 ) / ( r 2 − r 1 ) , aspect ratio of the stator
A R R = 1 2 ( H 3 + H 4 ) / ( r 4 − r 3 ) AR_{\mathrm{R}}=\tfrac{1}{2}(H_3 + H_4) / (r_4- r_3) A R R = 2 1 ( H 3 + H 4 ) / ( r 4 − r 3 ) , aspect ratio of the rotor
ξ S = ( h 2 − h 2 s ) / 0.5 v 2 2 \xi_\mathrm{S}= (h_2 - h_{2s})/0.5v_2^2 ξ S = ( h 2 − h 2 s ) /0.5 v 2 2 , loss coefficient of the stator
ξ R = ( h 4 − h 4 s ) / 0.5 w 4 2 \xi_\mathrm{R}= (h_4 - h_{4s})/0.5w_4^2 ξ R = ( h 4 − h 4 s ) /0.5 w 4 2 , loss coefficient of the rotor
Z S = 0.8 Z_\mathrm{S}= 0.8 Z S = 0.8 , Zweiffel parameter of the stator
Z R = 0.8 Z_\mathrm{R}= 0.8 Z R = 0.8 , Zweiffel parameter of the rotor
Outputs
The flow coefficient ϕ = v m / u 3 \phi=v_{m}/u_3 ϕ = v m / u 3 is computed from the stage equations
The flow angles β 3 \beta_3 β 3 and β 4 \beta_4 β 4 are computed from the stage equations
The spouting velocity is given by
v 0 = 2 ( h 01 − h 4 s ) v_0 = \sqrt{2(h_{01} - h_{4s})}
v 0 = 2 ( h 01 − h 4 s ) The blade velocity at each station is given by
u 1 = 0 u 2 = 0 u 3 = ( r 3 / r 4 ) u 4 = ν v 0 \begin{gather}
u_1 = 0 \\
u_2 = 0 \\
u_3 = (r_3/r_4) \\
u_4 = \nu \, v_0
\end{gather}
u 1 = 0 u 2 = 0 u 3 = ( r 3 / r 4 ) u 4 = ν v 0 The constant meridional velocity is given by:
v m = ϕ u 4 v_m = \phi \, u_4
v m = ϕ u 4 The absolute velocities are given by:
v 1 = v m / cos ( α 1 ) v 2 = v m / cos ( α 2 ) v 3 = v m / cos ( α 3 ) v 4 = v m / cos ( α 4 ) \begin{gather}
v_1 = v_m / \cos(\alpha_1) \\
v_2 = v_m / \cos(\alpha_2) \\
v_3 = v_m / \cos(\alpha_3) \\
v_4 = v_m / \cos(\alpha_4) \\
\end{gather}
v 1 = v m / cos ( α 1 ) v 2 = v m / cos ( α 2 ) v 3 = v m / cos ( α 3 ) v 4 = v m / cos ( α 4 ) The relative velocities are given by:
w 1 = v m / cos ( β 1 ) w 2 = v m / cos ( β 2 ) w 3 = v m / cos ( β 3 ) w 4 = v m / cos ( β 4 ) \begin{gather}
w_1 = v_m / \cos(\beta_1) \\
w_2 = v_m / \cos(\beta_2) \\
w_3 = v_m / \cos(\beta_3) \\
w_4 = v_m / \cos(\beta_4) \\
\end{gather}
w 1 = v m / cos ( β 1 ) w 2 = v m / cos ( β 2 ) w 3 = v m / cos ( β 3 ) w 4 = v m / cos ( β 4 ) The enthalpy of each station is given as
h 1 = h 01 − v 1 2 / 2 h_1 = h_{01} - v_{1}^2/2
h 1 = h 01 − v 1 2 /2 The enthalpy at the exit of the stator os obtained from the equation for the conservation of energy:
h 1 − h 2 u 3 2 = 1 2 ϕ 2 ( tan 2 α 2 − tan 2 α 1 ) \frac{h_1 - h_2}{u_3^2} = \frac{1}{2} \phi^2 \left( \tan^2 \alpha_2 - \tan^2 \alpha_1 \right)
u 3 2 h 1 − h 2 = 2 1 ϕ 2 ( tan 2 α 2 − tan 2 α 1 ) In the absence of losses and for constant meridional velocity, the enthalpy does not change across the interspace:
h 3 = h 2 , α 3 = α 2 h_3 = h_2, \quad \alpha_3 = \alpha_2
h 3 = h 2 , α 3 = α 2 The enthalpy at the exit of the rotor is obtained from the conservation of rothalpy:
h 3 − h 4 u 3 2 = 1 2 ϕ 2 ( tan 2 β 4 − tan 2 β 3 ) − 1 2 ( 1 − ρ 2 ) \frac{h_3 - h_4}{u_3^2} = \frac{1}{2} \phi^2 \left( \tan^2 \beta_4 - \tan^2 \beta_3 \right)- \frac{1}{2}(1-\rho^2)
u 3 2 h 3 − h 4 = 2 1 ϕ 2 ( tan 2 β 4 − tan 2 β 3 ) − 2 1 ( 1 − ρ 2 ) Once the enthalpy at each flow station is known, the thermodynamic state can be calculated assuming an isentropic process accross the turbine
[ p 1 , ρ 1 , a 1 , q 1 ] = E O S ( h 1 , s 1 ) [ p 2 , ρ 2 , a 2 , q 2 ] = E O S ( h 2 , s 2 ) [ p 3 , ρ 3 , a 3 , q 3 ] = E O S ( h 3 , s 3 ) [ p 4 , ρ 4 , a 4 , q 4 ] = E O S ( h 4 , s 4 ) \begin{gather}
[p_1, \rho_1, a_1, q_1] = \mathrm{EOS}(h_1, s_1) \\
[p_2, \rho_2, a_2, q_2] = \mathrm{EOS}(h_2, s_2) \\
[p_3, \rho_3, a_3, q_3] = \mathrm{EOS}(h_3, s_3) \\
[p_4, \rho_4, a_4, q_4] = \mathrm{EOS}(h_4, s_4) \\
\end{gather}
[ p 1 , ρ 1 , a 1 , q 1 ] = EOS ( h 1 , s 1 ) [ p 2 , ρ 2 , a 2 , q 2 ] = EOS ( h 2 , s 2 ) [ p 3 , ρ 3 , a 3 , q 3 ] = EOS ( h 3 , s 3 ) [ p 4 , ρ 4 , a 4 , q 4 ] = EOS ( h 4 , s 4 ) The equations for conservation of mass are used in order to compute all the radii and blade heights at all stations. From the mass balance in the stator we find that:
H 2 H 1 = ρ 1 ρ 2 r 1 r 2 \frac{H_2}{H_1} = \frac{\rho_1}{\rho_2} \, \frac{r_1}{r_2}
H 1 H 2 = ρ 2 ρ 1 r 2 r 1 from the definition of aspect ratio of the stator
H 2 r 2 = 2 A R S ( 1 − r 1 / r 2 ) ( 1 + H 1 / H 2 ) − 1 \frac{H_2}{r_2} = 2 AR_{\mathrm{S}}( 1 - r_1/r_2)(1 + H_1/H_2)^{-1}
r 2 H 2 = 2 A R S ( 1 − r 1 / r 2 ) ( 1 + H 1 / H 2 ) − 1 Inserting these results into the definition of mass flow rate we find that
m ˙ = ρ 2 v m A 2 = ρ 2 v m 2 π r 2 H 2 m ˙ = ρ 2 v 2 2 π r 2 2 [ A R S ( 1 − r 1 / r 2 ) ( 1 + H 1 / H 2 ) − 1 ] r 2 2 = m ˙ 4 π v m ρ 2 A R S − 1 ( 1 − r 1 / r 2 ) − 1 ( 1 + H 1 / H 2 ) \begin{gather}
\dot{m} = \rho_2 v_m A_2 = \rho_2 v_m 2\pi r_2 H_2 \\
\dot{m} = \rho_2 v_2 2\pi r_2^2 \left[AR_{\mathrm{S}} \, ( 1 - r_1/r_2)(1 + H_1/H_2)^{-1} \right] \\
r_2^2 = \frac{\dot{m}}{4\pi v_m \rho_2} AR_{\mathrm{S}}^{-1} \, ( 1 - r_1/r_2)^{-1}(1 + H_1/H_2)
\end{gather}
m ˙ = ρ 2 v m A 2 = ρ 2 v m 2 π r 2 H 2 m ˙ = ρ 2 v 2 2 π r 2 2 [ A R S ( 1 − r 1 / r 2 ) ( 1 + H 1 / H 2 ) − 1 ] r 2 2 = 4 π v m ρ 2 m ˙ A R S − 1 ( 1 − r 1 / r 2 ) − 1 ( 1 + H 1 / H 2 ) Once r 2 r_2 r 2 is known, the radius at the other stations is given by:
r 1 = r 2 ⋅ ( r 1 / r 2 ) r 3 = r 2 ⋅ ( r 2 / r 3 ) − 1 r 4 = r 3 ⋅ ( r 3 / r 4 ) − 1 \begin{gather}
r_1 = r_2 \cdot (r_1 / r_2) \\
r_3 = r_2 \cdot (r_2 / r_3)^{-1} \\
r_4 = r_3 \cdot (r_3 / r_4)^{-1}
\end{gather}
r 1 = r 2 ⋅ ( r 1 / r 2 ) r 3 = r 2 ⋅ ( r 2 / r 3 ) − 1 r 4 = r 3 ⋅ ( r 3 / r 4 ) − 1 The blade height at each station is thus given by:
H 1 = m ˙ ( 2 π r 1 ρ 1 v m ) − 1 H 2 = m ˙ ( 2 π r 2 ρ 2 v m ) − 1 H 3 = m ˙ ( 2 π r 3 ρ 3 v m ) − 1 H 4 = m ˙ ( 2 π r 4 ρ 4 v m ) − 1 \begin{gather}
H_1 = \dot{m} (2\pi r_1 \rho_1 v_m)^{-1} \\
H_2 = \dot{m} (2\pi r_2 \rho_2 v_m)^{-1} \\
H_3 = \dot{m} (2\pi r_3 \rho_3 v_m)^{-1} \\
H_4 = \dot{m} (2\pi r_4 \rho_4 v_m)^{-1} \\
\end{gather}
H 1 = m ˙ ( 2 π r 1 ρ 1 v m ) − 1 H 2 = m ˙ ( 2 π r 2 ρ 2 v m ) − 1 H 3 = m ˙ ( 2 π r 3 ρ 3 v m ) − 1 H 4 = m ˙ ( 2 π r 4 ρ 4 v m ) − 1 In order to compute the number of blades and openng of the stator blades we can use the Zweiffel criterion to determine the spacing-to-chord ratio
s S r 2 − r 1 = 0.5 Z S cos 2 ( α 2 ) ( tan α 1 − tan α 2 ) N S = π ( r 1 + r 2 ) / s S o S = s S cos ( α 2 + 1 2 2 π N S ) \begin{gather}
\frac{s_{\mathrm{S}}}{r_2- r_1} = \frac{0.5 Z_{\mathrm{S}}}{\cos^{2}\!\left( \alpha_{\mathrm{2}} \right) \left( \tan \alpha_{\mathrm{1}} - \tan \alpha_{\mathrm{2}} \right)} \\
N_{\mathrm{S}} = \pi (r_1 + r_2) / s_{\mathrm{S}} \\
o_{\mathrm{S}} = s_{\mathrm{S}} \cos(\alpha_2 + \tfrac{1}{2}\tfrac{2\pi}{N_{\mathrm{S}}})
\end{gather}
r 2 − r 1 s S = cos 2 ( α 2 ) ( tan α 1 − tan α 2 ) 0.5 Z S N S = π ( r 1 + r 2 ) / s S o S = s S cos ( α 2 + 2 1 N S 2 π ) Similarly, for the rotor we have that:
s R r 4 − r 3 = 0.5 Z R cos 2 ( β 4 ) ( tan β 3 − tan β 4 ) N R = π ( r 4 + r 3 ) / s R o R = s R cos ( β 4 + 1 2 2 π N R ) \begin{gather}
\frac{s_{\mathrm{R}}}{r_4- r_3} = \frac{0.5 Z_{\mathrm{R}}}{\cos^{2}\!\left( \beta_{\mathrm{4}} \right) \left( \tan \beta_{\mathrm{3}} - \tan \beta_{\mathrm{4}} \right)} \\
N_{\mathrm{R}} = \pi (r_4 + r_3) / s_{\mathrm{R}} \\
o_{\mathrm{R}} = s_{\mathrm{R}} \cos(\beta_4 + \tfrac{1}{2}\tfrac{2\pi}{N_{\mathrm{R}}})
\end{gather}
r 4 − r 3 s R = cos 2 ( β 4 ) ( tan β 3 − tan β 4 ) 0.5 Z R N R = π ( r 4 + r 3 ) / s R o R = s R cos ( β 4 + 2 1 N R 2 π )